Selective transmission of Dirac electrons and ballistic magnetoresistance of n-p 

junctions in graphene. 
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We show that an electrostatically created n-p junction separating the electron and hole gas regions 
in a graphene monolayer transmits only those quasiparticles that approach it almost perpendicularly 
to the n-p interface. Such a selective transmission of carriers by a single n-p junction would manifest 
itself in non-local magnetoresistance effect in arrays of such junctions and determines the unusual 
Fano factor in the current noise universal for the n-p junctions in graphene. 
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The chiral nature of quasiparticles in graphene mono- 
layers and bilayers (i, iS, ^ li BL has been revealed in 
several recent experiments [3,11, E, 0|- The Fermi 
level in a neutral graphene sheet (a monolayer of carbon 
atoms with hexagonal lattice structure) is pinned near 
the corners of the hexagonal Brillouin zone which 
determine two non-equivalent valleys in the quasiparticle 
spectrum. The quasiparticles in each of the two valleys 
are described by the Hamiltonian [^Q, 

Hi = V a ■ p, 

where the 'isospin' Pauli matrices ai operate in the space 
of the electron amplitude on two sites (A and B) in the 
unit cell of a hexagonal crystal p = (px^Py) — —iV 
is the momentum operator [T^ defined with respect to 
the centre of the corresponding valley, and w is a con- 
stant formed by the A-B hopping 4]. The Dirac-type 
Hamiltonian Hi determines the linear dispersion vp for 
the electron, and —vp for the hole branch of quasiparti- 
cles. In each valley jlJi, the electron and hole states also 
differ by the isospin projection onto the direction of their 
momentum: electrons have chirality a ■ p/p = 1, holes 

■ p/p — ^1- Therefore, in structures where the quasi- 
particle isospin is conserved (, a monolayer with electro- 
static potential scattering) their backscattering is strictly 
forbidden |^, which gives rise to the pecuhar properties 
of the n-p junction in graphene reported in this Letter. 

Since an atomically-thin graphitic film is a gapless 
semiconductor, carrier density in it can be varied u sing 
external gates ^ from electrons to holes 0, H, 0, E ^M- 
A planar n-p junction in graphene can be made, e.g., 
using split-gates, and in view of a rapidly improving mo- 
bility of the new material d 0, it may soon be possi- 
ble to fabricate ballistic circuits of electrically controlled 
graphene-based n-p junctions. Below, we model the n- 
p junction in graphene using the electrostatic potential 
u{x) = vkFr](x / d) characterized by a single length scale 
d and the Fermi momentum kp determined by the equal 
densities of the electron and hole gases on the opposite 
sides of it. Here 77(±oo) — ±1, ?7'(0) — 1, and the line 
X — separates the p- and n-regions. Since in a junction 
produced by electrostatic gates the length d is about the 
inter-gate distance and exceeds the electron wavelength 
in a monolayer, we focus this study on smooth n-p junc- 



tions with kpd > 1, and show that their transmission 
properties are determined by the central region where 
u{x) ^ Fx[F = vkp/d]. 

The transport properties of the n-p junction are de- 
termined by the angular dependence of the probability 
w{9) of an electron incident from the left with an en- 
ergy equal to the chemical potential fi = to emerge 
as a hole on the right hand side of the junction. The 
electron approaching the center of the junction has the 

kinetic energy 



-p"^, where p„ = kf 



Py = kp sm a IS con- 
served. Due to the energy conservation, the x-component 
of the electron momentum is Px{x) = ^u'^{x)/v'^ — Py, 

the classically allowed region for the electron motion is 
determined by the condition \u\ > pyV and its trajectory 
cannot extend beyond the turning point at the distance 
I = vpy/F from the center of the junction For a 

particle incident perpendicular to the junction (py — 0) 
the classically forbidden region disappears. Moreover, 
due to the isospin conservation which prohibits backscat- 
tering of chiral quasiparticles J^, the wave incident at 
= is perfectly transmitted, though, for any small 9, 
the transmission probability is determined by tunnelling 
through the classically forbidden region, w ~ e~^'^, 
where S — i J ^px{x)dx — ^irvpy/F. For a smooth n-p 
junction shown in Fig. ^with F = vkp/d and kpd 3> 1, 
this yields (for 6* < 7r/2) 



w{e) 



-7r(fci?rf) sin 9 



(1) 



The angular dependence of the transmission probabil- 
ity given in Eq. ^ is, in fact, exact in the range 9 ^ 7r/2 
for any smooth junction and represents the central result 
of this Letter. Below, we rigorously derive the result in 
Eq. Q using the method of transfer matrix and also 
show that w{9) = cos^ 9 for a step-like potential. Sim- 
ilarly to the formulae describing adiabatic ballistic 
constrictions in semiconductors, the applicability of Eq. 

is not restricted by the constraint w <C 1. This can be 
used to describe how a smooth n-p junction selectively 
transmits only carriers approaching it within a small an- 
gle 9 < 00 — l/y/nkpd around the perpendicular direc- 
tion and to determine the conductance per unit length of 
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FIG. 1: Angular dependence of quasiparticle transmis- 
sion through the electrostatically generated n-p junction in 
graphene. 



a broad junction, 



9np- ^ 



kpd9 2e^ kp 

2tt nh V d ' 



and the universal Fano factor in the shot noise, 



((/•/)) = (l-./i)e/ 



(2) 



(3) 



At the end of this paper we shall discuss several ballistic 
magnetoresistance effects which exploit the selectivity of 
transmission implicit in Eq. 

To formulate the scattering problem, we shall exploit 
the separation of x and y variables for the electron mo- 
tion across the junction (in x direction) and the fact 
that momentum along y-axis (parallel to the junction) 
is conserved. This makes the scattering problem one- 
dimensional (ID). The scattering states at the energy 
equal to the chemical potential, /i = are spinors satis- 
fying the Dirac-type equation 



-idxTx^ + V ^u{x)^ + Pyaytp — 0, 



(4) 



which conserves the ID current = ip'^cTxtp. 

To find the transmission probability w{9) for such 
states, we calculate the transfer matrix T{x, y) which sat- 
isfies the equation 



dxT{x,y) = L{x)T{x,y), L 



.u{x) 

-I (Tt 



+ Py(yz (5) 



and the conditions T{y,y) — I, T{x,y) — T{x, z)T{z,y), 
det T{x,y) = 1, and T\x,y)axT{x,y) = a^- To re- 
late the transmission coefficient w to the transfer matrix 
T{x,y) one has to factor out the asymptotic evolution of 
the reflected and transmitted waves. This can be done 
by using matrices A± satisfying the wave equation in the 
asymptotic regions, 



dxA±{x) = {Tikpo-x +Py<Jz)A±{x), 



(6) 



such that their columns are made of right- and left- 
propagating states normalized to carry unit current. The 
explicit expression for these matrices is 



A±{x) 



2^ V e'F'^^"'' 



r,±ip^X 



where Px{x) = yu'^{x)/v^ ^ Py ~ 
transmission probability can be found using the matrix 



k'p 



Py. Then, the 



a (3* 



lim A-^{x)T{x,-x)A_{-x), 



1 



w 



(7) 

To illustrate the transfer matrix formalism, we cal- 
culate the probability of a Dirac fermion transmission 
through a sharp potential step u{x) = vkpsigriix). In 
this case, we factor the transfer matrix as T(x^ y) = 
r+(a;, 0)T_(0, y), where T+(_)(a;,?/) is a transfer matrix 
on the right(left) side of the junction, each given by 
T±{x,y) — A±{x)A'^^ {y). Using this solution and Eq. 
O we find 



a P* 
(3 a * 



= A;i(0)A_(0), a = l- 



iPy 



k^ 



For the transmission probability this yields 
■w = 1 — {py/kpY — cos^ 9, 



(8) 



which manifests the chiral nature of quasiparticles in 
graphene. Indeed, the free electron states of the 
Dirac Hamiltonian Hi have their isospin polarized along 
the momentum (for the transmitted holes, with p = 
{—kp cos9, kp sm9), it is antiparallel) , and the reflection 
amplitude of an electron is determined by the scalar prod- 
uct ip'y, ■ t/j'^ ^ sm9 of its initial and flnal state spinors. 

To calculate the transmission probability for a smooth 
potential with kpd ^ 1, we separate the x-axis across 
the junction into the inner (i) and outer (o) parts. In the 
outer part, \x\ > cd (where c <C 1), we flnd the T-matrix, 
To using the method of adiabatic expansion. Then, we 
match it with the exact solution, Ti obtained in the cen- 
tral part of the junction, |a:| < d, where the potential u{x) 
can be linearised, u{x) w kpx/d, and obtain the complete 
tranfer matrix as T{y,x) = To{x,a)Ti{a,—a)To{—a,y). 

For the adiabatic expansion of the transfer matrix To 
we use a transformation 



Y{x) 



IK IX 

u \ u u 



p-i"") ^ y ^ - pI (9) 

K^Py + ipx{x). 



which locally diagonaliscs the L operator in Eq. ([SJ): 

Y-^LY ^ipx{x)a,. (10) 
The transfer matrix To defined in a new basis, 

To{x,y)=Y{x)fo{x,y)Y-\y), (11) 
satisfies the equation 

dxfo{x,y) = ipxix)a^foix,y) + fl{x)fo{x,y), (12) 



= -Y-'dxY = 



Pyu'{x) f —H yc' 



2pl{x)u{x) 
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In the adiabatic approximation the matrix fl(x) is as- 
sumed to be small as compared to the diagonal term 
Px{x)<Jz, and to the leading order Eq. ((T^ is solved by 



fo{x,y) = exp 



icfz / Px{x')dx' 



(13) 



Formally, the adiabatic approximation is justified if 
\pyu' / {up'^)\ <C 1, which breaks down near the turning 
points Px{x) — and when u{x) = 0. However, for the 
junctions with kpd 3> 1, the interval between turning 
points lies within the region of space where the the poten- 
tial profile can be approximated using the linear function 
u{x) = kpx/d. The transfer matrix in this region, Ti can 
be found from Eq. Q exactly, using the transformation 



T,{x,y) = e-^^'^='e-'^'^^f,(x,y)e'^^»e'^'^«. 



(14) 



where (j){x) = kpd ^x^. 
This is because the matrix Ti satisfies the equation 



dxf^{x,y) = -py 







T^ix,y), (15) 



where the upper row of Ti can be expressed in terms of 
two linearly independent solutions of the equation 

while the lower row can be expressed in terms of their 
complex conjugate. Equation H15|l is symmetric with 
respect to the parity transformation x — > —x, and its 
even/odd solutions are 



^odd{x)^~PyX<i>(^^-i^, 



where $ is the confluent hypergeometric (Kummer) func- 
tion ^3 with the following asymptotic properties: 



$(a, b;z ^ ioo) 



1 



T{b-a) z"" r(a) 



Therefore, inside the interval |a;|,|y| < cd the transfer 
matrix T,- can be written as 



f,{x,y) = B{x)B~\y), B 



odd oven 



(16) 



where the matrix B satisfies Eq. (|15|) and has unit Wron- 
skian, det B = 1. 

Finally, after a chain of substitutions, the obtained so- 
lutions for the matching transfer matrices To and Ti can 
be combined together into 

T{y,x) = To{x,a)T^{a,-a)Toi-a,y), 



and used to calculate the parameters a and f3 in Eq. 10, 



a — e ^''F , 



r(i + l^) 

[px{x') - px{oo)]dx' , 

needed for determining the transmission probability, 

w = \a\^^ = e-''P>^^^. (17) 

A selective transmission of carriers by a smooth n- 
p junction described by Eqs. H17ll|l . with kpd ^ 1, 
only allows for the passage of quasiparticlcs approaching 
the junction in an almost perpendicular direction, with 
Py < ^Jkpjd <^ kp and 0^1. This makes the transport 
characteristics of ballistic graphene-based devices sensi- 
tive to the geometrical orientation of n-p junctions in 
them, and it is capable of generating a sizeable magne- 
toresistance (MR) effect. 

A nominal resistance, R^p — 1/agnp of a single, sep- 
arately taken n-p junction with the peripheral length 
a separating the electron and hole gases with densities 
"-e//i = kp/T^ is determined by Eq. 0. Whether or 
not the nominal junction resistance contributes to the 
total resistance of a ballistic device depends on how free 
carriers propagate in it. For example, when an n-p junc- 
tion, with the perimetre a = 27rr, separates two metallic 
Corbino contacts to the ballistic 2D electron/hole gases 
shown in Fig. |2Ia), electrons emitted from the inner 
contact with the radius b < r/\/nkpd reach the junction 
at the incidence angle 6 < 9q = l/^/nkpd and pass it 
without scattering. As a result, the presence of the n-p 
junction does not affect the Corbino resistance, unless 
an external magnetic field changes the incidence angle to 



0' = r/rc > 9o, where 



kphc/eB is the cyclotron 



radius in the ballistic region. This generates the MR, 
R.{B) = i?cxt + where /(O) - 0, /(I) - 1, 



and = {hc/e) ^Jhrj-KT'^d. 



(18) 



A strong MR effect can also be expected an a 
Hall-bar sample with several parallel n-p-n junctions. 
Fig. I21b). The energy-averaged [l^ transmis- 
sion through the series of two junctions, W2{0) = 

[w^^{9) + w^^{9 + r/rc) — l] is determined by the in- 
dividual junction transmissions 111(0) a,ndw{9-\--^). Here, 
we take into account that, due to the external magnetic 
field, an electron transmitted by the first junction at the 
incidence angle 9 would approach the second at the an- 
gle 6' ^ 9 + ^ , where r/rc ^ B/B^, with defined in 
Eq. (|18|l . In the absence of a field 9' = 9, and the trans- 
mitted particle would also pass the second junction, as 
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FIG. 2: Ballistic MR devices with n-p junctions in graphene: 
(a) Corbino geometry; (b) series of n-p-n junctions, with the 
illustration of trajectories of electrons transmitted by the first 
junction for B — (left) and B > B* (right); (c) three- 
terminal cavity. 

shown on the left of Fig. |2Ib). If, due to a high field, 
the angle 9' is sufficient for the particle to be refiected, 
9' > 9q the latter would return to the first junction along 
the path illustrated on the right in Fig. |2Ib) and escape 
to the contact where it came from. This would suppress 



the conductance of the n-p-n junction down to the value 
determined scattering by the side edges of the sample. 
Having substituted W2{9) [instead of w{9)] into the con- 
ductance per unit length of a broad junction defined in 
Eq. (0), we find the magnetoconductance of the n-p-n 
j miction, 

g„,„(i?) « ^ /°° , (ff... ■ (19) 

A strongly selective quasiparticle transmission in Eqs. 
(|1I17|) can also be used for creating ballistic cavity-type 
structures in graphene, with non-local transport prop- 
erties. In a three-terminal 'cavity' shown in Fig. EJc), 
a p-charging gate would produce two parallel n-p junc- 
tions, so that ballistic electrons emitted from the contact 
1 and transmitted by the first junction would easily pass 
through the second and reach contact 3. As a result, 
a bias voltage applied between contacts 1 and 2 would 
generate current between contacts 1 and 3, thus giving 
rise to the trans-conductance Gjf with a strong magnetic 
field dependence, 

asm^'-^^n^i (20) 
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